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Abstract 

A  parametric  analysis  of  a  model  equation  developed  to  predict  the  cell  voltage  versus  current  density  response  of  a  liquid  feed  direct  methanol 
fuel  cell  is  presented.  The  equation  is  based  on  a  semi-empirical  approach  in  which  methanol  oxidation  and  oxygen  reduction  kinetics  are  combined 
with  effective  mass  transport  coefficients  for  the  fuel  cell  electrodes.  The  model  equation  is  applied  to  experimental  data  for  a  small-scale  fuel 
cell  and  produces  electrochemical  parameters  generally  consistent  with  those  expected  for  the  individual  components  of  the  fuel  cell  MEA.  The 
parameters  thus  determined  are  also  used  in  the  model  to  predict  the  performance  of  a  DMFC  with  a  new  membrane  electrode  assembly. 
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1.  Introduction 

Previous  models  of  the  DMFC  have  been  few.  A  one¬ 
dimensional  model  of  the  potential  distribution  and  concentra¬ 
tion  distribution  of  methanol  in  the  anode  electrocatalyst  layer 
was  used  in  a  model  of  the  overall  cell  voltage  performance  of  a 
vapour  feed  system  [1,2].  The  model  gives  good  agreement  with 
experimental  data  except  under  conditions  where  mass  transport 
becomes  rate  limiting.  Baxter  et  al.  [3]  have  presented  a  model 
of  the  DMFC  anode  which  is  considered  to  be  a  porous  electrode 
consisting  of  an  electronically  conducting  catalyst  structure  that 
is  thinly  coated  with  an  ion- selective  polymer  electrolyte.  The 
pores  of  the  electrode  are  filled  with  aqueous  methanol  solu¬ 
tion  in  which  all  species  of  the  reaction  are  free  to  transport. 
Mass  transfer  in  the  anode  is  defined  in  terms  of  a  pseudo-mass 
transport  coefficient.  Other  models  have  ignored  current  and 
potential  distributions  occurring  in  the  electrocatalyst  structures 
and  considered  the  mass  transport  processes  of  methanol  in  the 
cell  and  the  effect  on  electrode  performance  [4].  Recently  we 
have  considered  the  applicability  of  empirical  models,  originally 
developed  for  hydrogen  polymer  electrolyte  fuel  cells  (PEMFC) 
[5]. 
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The  main  drawback  of  the  empirical  equations  derived  for 
the  PEMFC  [6-9]  was  their  inability  to  explain  theoretically  the 
values  of  parameters  used  to  fit  the  data.  In  the  first  part  of  this 
paper  [10]  we  confirmed  the  general  ability  of  the  semi-empirical 
model  to  predict  experimental  data.  In  this  paper  we  consider 
the  components  of  the  model  in  detail  and  perform  a  parametric 
analysis  of  the  coefficients  obtained.  The  model  is  then  used  to 
predict  the  behaviour  of  a  new  DMFC. 

2.  Experimental 

Tests  on  the  DMFC  were  performed  with  a  cell  with  a  cross- 
sectional  area  of  9  cm2.  The  cell  consisted  of  two  non-porous 
graphite  blocks  with  a  series  of  parallel  channels  machined  into 
the  block  for  the  flow  of  methanol  and  oxygen/air  (Fig.  1). 

The  cell  was  fitted  with  a  membrane  electrode  assembly 
(MEA)  sandwiched  between  the  two  graphite  blocks.  The  cell 
was  held  together  between  two  plastic  insulation  sheets  and 
two  stainless  steel  backing  plates  using  a  set  of  retaining  bolts 
positioned  around  the  periphery  of  the  cell.  Electrical  heaters, 
supplied  by  Watson  Marlow,  were  placed  behind  each  of  the 
graphite  blocks  to  heat  the  cell  to  the  desired  operating  tem¬ 
perature.  The  graphite  blocks  were  also  provided  with  electrical 
contacts  and  small  holes  to  accommodate  thermocouples.  The 
fuel  cells  were  used  in  a  simple  flow  rig,  which  consisted  of 
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Fig.  1.  DMFC  small-scale  cell  flow  bed  design  and  manifold  arrangement. 


a  Watson  Marlow  peristaltic  pump  to  supply  aqueous  methanol 
solution,  from  a  reservoir,  to  a  Eurotherm  temperature  controller 
to  maintain  the  cell  at  a  constant  temperature.  Air  was  supplied 
from  cylinders,  at  ambient  temperature,  and  the  pressure  regu¬ 
lated  by  needle  valves. 

MEAs  studied  in  this  work  were  made  in  the  following  man¬ 
ner:  the  anode  consisted  of  a  Teflonised  (20%)  carbon  cloth 
support  (E-Tek,  type  ‘A’),  of  0.3  mm  thickness,  upon  which 
was  spread  a  thin  (diffusion  layer)  layer  of  uncatalysed  (ket- 
jenblack  600)  10wt.%  Teflonised  carbon.  The  catalysed  layer, 
unless  otherwise  specified  consisted  of  35  wt.%  Pt-15  wt.%  Ru 
(2  mg  cm-2  metal  loading)  dispersed  on  carbon  (ketjen)  and 
bound  with  bound  with  10  wt.%  Nation®  from  a  solution  of 
5  wt.%  Nation®  dissolved  in  a  mixture  of  water  and  lower 
aliphatic  alcohol’s  (Aldrich),  was  spread  on  this  diffusion  back¬ 
ing  layer.  A  thin  layer  of  Nation®  solution  was  spread  onto  the 
surface  of  each  electrode.  The  cathode  was  constructed  using 
a  similar  method  as  for  the  anode,  using  a  thin  diffusion  layer 
bound  with  10  wt.%  PTFE,  and  1  mg  cm-2  Pt  black  with  10  wt.% 


Nation®  in  the  catalyst  layer.  The  electrodes  were  placed  either 
side  of  a  pre-treated  Nation®  117  membrane.  This  pre-treatment 
involved  boiling  the  membrane  for  1  h  in  5  vol.%  H2O2  and  1  h 
in  1  mol  dm-3  H2SO4  before  washing  in  boiling  Millipore  water 
(>18m£2)  for  2h  with  regular  changes  of  water.  The  assembly 
was  hot-pressed  at  100  kg  cm-2  for  3  min  at  135  °C. 

Cell  voltage  versus  current  density  response  was  measured 
galvanostatically,  by  incrementally  increasing  the  current  from 
open  circuit  and  measuring  the  cell  potential  and  then  reducing 
the  current  incrementally  again  measuring  the  cell  voltage.  The 
data  reported  here  were  obtained  with  one  MEA,  unless  other¬ 
wise  stated.  The  MEA  was  conditioned  before  use  in  two  stages: 
for  48  h  in  the  test  cell,  in  2.0  mol  dm-3  methanol  solution  at 
75  °C,  and  then  by  maintaining  the  cell  with  an  applied  load  of 
100  mA  cm-2  for  several  hours.  This  pre-treatment  resulted  in 
stable  performance  under  continuous  operation. 

3.  Model  equation 


As  a  kinetic  model  for  methanol  oxidation,  Tafel  type  kinetics 
is  chosen: 


j  =  jo 
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where  jo  is  the  exchange  current  at  the  reference  concentration, 
a  the  transfer  coefficient,  “a”  refers  to  the  condition  at  the  anode 
catalyst  surface  and  N  is  an  order  of  reaction. 

We  represent  mass  transport  in  the  anode  side  of  the  cell  using 
the  effective  mass  transport  coefficient,  &eff- 


j  =  keffnF(CME  ~  C^e) 
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By  re-arranging  the  previous  Eq.  (1)  the  overpotential 
becomes: 
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By  combining  Eqs.  (2)  and  (3)  we  obtain: 
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Using  Eq.  (4)  for  the  anode  we  can  obtain  an  expression  for  the 
cell  voltage  by  introducing  the  kinetic  expression  for  oxygen 
reduction: 


j  =  70c 
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where  po  is  the  partial  pressure  for  oxygen  and  No  is  the  reaction 
order  for  oxygen  reduction. 

The  cathode  overpotential  is  then 
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where  k\o  is  the  mass  transport  coefficient  for  the  cathode  side 
of  the  cell. 

Combining  Eqs.  (4)  and  (6),  we  obtain  an  expression  for  the 
cell  voltage,  Ece\\,  as 
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To  model  the  performance  of  the  DMFC  we  propose  the  fol¬ 
lowing  semi-empirical  equation  for  the  cell  voltage,  £ceii,  in 
which  we  assume  that  the  reduction  of  oxygen  does  not  proceed 
under  mass  transport  limitations: 

£ceii  =  E*o  -  /’cell  log  j  -  Rej  +  Cl  ln(l  -  C2J)  (8) 


Cell  Temperature/  K 


Fig.  2.  Fitted  values  for  parameter  E*0  of  Eq.  (8)  as  a  function  of  cell  temperature 
in  K  for  a  cell  operated  with  various  methanol  solutions  (methanol  solution 
concentrations,  A:  0. 125  M;  0:  0.25  M;  x :  0.5  M;  0:  0.75  M)  supplied  at  a  rate 
of  1.12  cm3  min-1  with  air  fed  cathodes  pressurised  at  2  bar. 
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In  the  above  j  is  the  current  density,  aa  and  ac  the  transfer  coef¬ 
ficients  for  the  oxidation  of  methanol  and  reduction  of  oxygen, 
respectively,  R&  the  internal  cell  resistance  (mainly  due  to  the 
polymer  electrolyte  membrane),  N  a  reaction  order  for  methanol 
oxidation,  &eff  an  effective  mass  transport  coefficient  for  the 
anode  side  of  the  cell,  Cme  the  methanol  concentration  and 
po  is  the  oxygen  pressure.  Cocell  is  the  standard  potential  for 
the  DMFC  overall  reaction,  which  theoretically  is  given  by  the 
Nemst  equation. 


4.  Application  of  the  model 


model  Eq.  (1).  Open  circuit  potential  also  increases  with  increas¬ 
ing  temperature  as  expected  due  to  more  catalyst  activity  and  as 
predicted  from  the  Nemst  equation.  These  are  general  trends  and 
a  more  detailed  theoretical  approach  can  be  found  below  in  the 
relative  section. 

Fig.  3  shows  the  combined  effect  of  methanol  solution  con¬ 
centration  and  temperature  on  the  effective  Tafel  slope  Z?ceii.  As  it 
can  be  seen  Tafel  slope  increases  with  an  increase  in  cell  temper¬ 
ature  and  is  not  influenced  by  methanol  solution  concentration  to 
a  significant  extent.  The  effect  of  temperature  is  consistent  with 
that  predicted  by  the  model.  For  example  an  increase  in  tem¬ 
perature  from  303  to  363  K  gives  an  approximate  20%  increase 
in  the  value  of  bce\\  (0.108-0.124).  It  should  be  noted  here  that 
an  increase  in  methanol  solution  concentration  is  not  expected 
to  effect  of  value  of  bCQ n  unless  of  course  then  is  an  influence 
brought  about  by  methanol  crossover. 

Fig.  4  shows  the  combined  effect  of  methanol  solution  con¬ 
centration  and  temperature  on  cell  resistance  R.  The  effect  of 


In  the  previous  paper  [10]  we  confirmed  the  general  applica¬ 
bility  of  the  model  Eq.  (1)  to  represent  the  cell  voltage,  current 
density  behaviour  of  the  DMFC. 

The  proposed  model  was  used  to  predict  the  voltage  response 
of  a  specific  small-scale  DMFC  cell  when  operated  under  con¬ 
stant  cathodic  overpressure  with  air  fed  at  a  constant  rate, 
and  with  the  anode  side  solution  flow  rate  also  constant  at 
1.12  cm3  min-1.  In  doing  so  we  tried  to  minimise  the  effects 
that  are  attributable  to  changes  in  the  cathode  side  overpressure 
or  due  to  an  increased  methanol  solution  flow  rate.  Both  param¬ 
eters  under  specific  circumstances  can  alter  the  cell  response 
[11-17]. 

Fig.  2  indicates  that  the  open  cell  voltage  decreases,  by 
approximately  0.05  V,  with  increasing  methanol  solution  con¬ 
centration  from  0.25  to  0.5  M,  as  the  impact  of  methanol 
crossover  on  cathode  performance  becomes  more  significant. 
The  variation  in  Eq  is  consistent  with  that  as  proposed  in  the 


Fig.  3.  Fitted  values  for  parameter  bce\\  of  Eq.  (1)  as  a  function  of  cell  temperature 
in  K  for  a  cell  operated  with  various  methanol  solution  concentrations,  supplied 
at  a  rate  of  1 . 12  cm3  min-1  with  air  fed  cathodes  pressurised  at  2  bar  (A:  0. 125  M; 
O:  0.25  M;  x:  0.5  M;  0:  0.75  M). 
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Cell  Temperature/  K 

Fig.  4.  Fitted  values  for  parameter  R  of  Eq.  (1)  as  a  function  of  cell  temperature 
in  K  for  a  cell  operated  with  various  methanol  solutions  (methanol  solution 
concentrations,  A:  0.125  M;  0:  0-25  M;  x :  0.5  M;  0:  0.75  M)  supplied  at  a  rate 
of  1.12  cm3  min-1  with  air  fed  cathodes  pressurised  at  2  bar. 

methanol  solution  concentration  on  resistance  R  is  negligible  in 
the  concentration  range  considered.  The  effect  of  temperature  is 
to  decrease  the  resistance  as  expected  from  known  correlations 
for  the  resistivity  of  Nation™  membranes.  The  resistivity  of 
Nation  over  the  temperature  range  considered  in  Fig.  4  is  con¬ 
sistent  with  published  data.  For  example  at  80  °C  the  protonic 
conductivity  is  0.202  S  cm-1  which  with  a  170  p,m  thick  mem¬ 
brane  gives  a  resistivity  of  approximately  0.085  Qcm2,  which 
compares  favourably  with  the  data  in  Fig.  4. 

The  values  of  resistivity  are  correlated  by  the  expression: 

(B  B\ 

*  =  *oexP  (9) 

Ro  =  0.085  Q  cm2  is  the  resistance  at  a  temperature 
To  =  383.15  K  and  #  =  3724  is  a  constant  determined  from 
the  existing  data  set  (see  dotted  line  in  Fig.  3).  Hence  in  the 
model  Eq.  (9)  the  resistance  parameter  can  be  fixed,  with 
a  function  following  the  trend  indicated  in  Fig.  3,  and  thus 
reduces  the  number  of  fitting  parameters  in  the  model.  It  should 
be  noted  that  part  of  the  cell  resistance  is  not  due  to  the  Nation 
membrane. 

Fig.  5  shows  the  combined  effect  of  methanol  solution  con¬ 
centration  and  temperature  on  the  model  parameter  C\ .  As  pre¬ 
dicted,  the  value  of  C\  increases  with  increasing  temperature. 
The  effect  of  methanol  concentration  is  to  produce  a  mini¬ 
mum  value  of  coefficient  C\  at  a  particular  temperature.  This 
behaviour  is  oxidation  is  not  a  simple  electrochemical  process 
and  that  the  Tafel  kinetic  equation  is  indeed  an  approximation.  It 
has  been  reported  that  the  reaction  order  for  methanol  is  difficult 
to  explain  in  terms  of  the  model  equation.  A  contributing  fac¬ 
tor  will  be  that  the  reaction  order  for  methanol  oxidation  is  1 .0 
at  low  concentrations  and  decreases  to  zero  at  higher  concen¬ 
trations  of  around  2.0  M.  The  variation  in  reaction  order  itself 
depends  upon  the  electrode  potential.  In  terms  of  actual  coef¬ 
ficient  values,  it  is  expected  that  these  would  be  closely  linked 
to  the  Tafel  slope  for  the  reaction  ( RT/aaF ).  For  example  at  a 
temperature  of  343  K  and  with  a  transfer  coefficient,  of  0.5  and  a 


Cell  Temperature/  K 


Fig.  5.  Fitted  values  for  parameter  C\  of  Eq.  (8)  as  a  function  of  cell  temperature 
in  K  for  a  cell  operated  with  various  methanol  solutions  (methanol  solution 
concentrations,  A:  0.125  M;  0:  0-25  M;  x :  0.5  M;  0:  0.75  M)  supplied  at  a  rate 
of  1.12  cm3  min-1  with  air  fed  cathodes  pressurised  at  2  bar. 
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Methanol  Solution  Concentration/  M 

Fig.  6.  Fitted  values  for  parameter  C2  of  Eq.  (1)  as  a  function  of  aqueous 
methanol  solution  concentrations  in  M  for  five  different  cell  temperatures. 

reaction  order,  N=  1.0,  the  coefficient  C\  has  a  value  of  approx¬ 
imately  0.06  V.  This  value  agrees  with  the  computed  (“best  fit”) 
value  obtained  from  the  data  analysis  as  shown  in  Section  6  later. 

Fig.  6  shows  the  combined  effect  of  methanol  solution  con¬ 
centration  and  temperature  on  model  parameter  C2.  From  the 
definition  of  this  parameter  it  is  expected  that  an  increase  in 
methanol  concentration  would  decrease  the  value  of  C2,  which 
is  confirmed  in  Fig.  5.  In  addition  the  observed  decrease  in  the 
value  of  C2  with  temperature  is  expected  due  to  an  increase  in 
the  value  of  mass  transport  coefficient  at  higher  temperature:  a 
diffusion-related  effect.  We  explore  the  relationship  between  the 
model  coefficients  and  kinetic  and  mass  transport  parameters  in 
more  detail  below. 

5.  Effective  mass  transfer  coefficients 

The  variation  in  effective  mass  transfer  coefficient  with  tem¬ 
perature  is  partly  due  to  an  increase  in  methanol  diffusion 
coefficient  with  temperature.  An  additional  effect  created  by 
the  higher  temperatures  may  be  enhanced  mass  transport  due 
to  bubble  liberation  at  the  surface  of  the  ME  A.  At  higher  tem¬ 
peratures  a  higher  bubble  velocity  in  the  diffusion  layers  can 
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Fig.  7.  Empirical  equation  and  proposed  model  predicted  effective  mass  transfer 
coefficients  as  a  function  of  cell  temperature  in  for  a  cell  operated  with  various 
methanol  solutions  (methanol  solution  concentrations,  •:  0.125  M;  ♦:  0.25  M; 
■:  0.5  M;  A:0.75  M)  supplied  at  a  rate  of  1.12  cm3  min-1  with  air  fed  cathodes 
pressurised  at  2  bar.  Model  predictions  are  presented  with  dotted  lines. 

lead  to  an  increase  in  liquid  voidage,  which  in  turn  increases 
the  effective  mass  transfer  in  the  carbon  cloth.  It  is  not  known 
whether  a  higher  rate  of  CO2  generation  (higher  current  density) 
increases  the  gas  velocity  in  the  cloth  or  increases  the  gas  voidage 
in  the  cloth.  In  addition  higher  temperatures  will  increase  the 
methanol  and  water  transfer  through  the  Nation  membrane, 
which  will  influence  the  mass  transport  rate  of  methanol  at  the 
anode. 

The  methanol  concentration  influences  mass  transport  in  sev¬ 
eral  ways;  at  low  methanol  concentrations,  limiting  currents  are 
lower  and  therefore  lower  gas  evolution  rates  result.  Lower  gas 
evolution  rates  will  reduce  any  hydrodynamic  influence  on  the 
mass  transport  rate  at  the  carbon  cloth,  or  locally  at  the  surface 
of  the  anode  catalyst.  The  gas  evolution  rate  will  influence  the 
voidage  of  the  liquid  in  the  carbon  cloth,  which  will  in  turn  affect 
the  methanol  diffusion  rate.  Lower  methanol  concentrations 
also  reduce  methanol  transfer  by  diffusion  and  electro-osmosis, 
through  the  membrane,  due  to  a  lower  current  density  of  oper¬ 
ation,  which  will  also  reduce  the  methanol  transfer  rate  to  the 
anode.  The  local  temperature  and  pressure  at  the  anode  will 
also  influence  mass  transport  rates.  Higher  temperatures  increase 
the  local  vaporisation  rate  of  methanol  into  carbon  dioxide  gas 
formed  at  the  anode,  which  will  reduce  the  effective  liquid  phase 
concentration. 

Fig.  7  shows  the  variation  of  effective  mass  transfer  coeffi¬ 
cient  with  cell  temperature  and  methanol  solution  concentration. 
Mass  transfer  coefficients  increase  with  an  increase  in  tem¬ 
perature  and  with  methanol  concentration.  The  latter  effect  is 
essentially  due  to  the  higher  limiting  current  densities  at  the 
higher  concentrations,  which  cause  a  greater  evolution  of  car¬ 
bon  dioxide  gas.  The  data  in  Fig.  7  can  be  correlated  by  the 
expression: 
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3.342514414 

A 

4601.967164 

4014.656868 

3675.977493 

4505.215133 

The  prediction  of  Eq.  (10)  are  compared  with  the  experimental 
variation  of  mass  transfer  coefficient  with  methanol  concentra¬ 
tion  as  a  function  of  cell  temperature,  in  Fig.  7. 

6.  Empirical  equation-based  kinetic  parameter 

Fig.  8  shows  the  variation  in  the  combined  transfer  coef¬ 
ficients  alpha,  aaac/(aa  +  ac),  as  obtained  from  the  value  of 
the  coefficient  Z?ceii.  As  can  be  seen  the  variation  in  the  coef¬ 
ficient  with  methanol  concentration  is  small.  The  value  indi¬ 
cates  that  the  transfer  coefficients  for  one  or  both  cell  reactions 
are  greater  than  0.5  and  potentially  close  to  1.0.  The  effect 
of  temperature  is  to  increase  the  value  of  alpha  by  approxi¬ 
mately  5%  over  the  range  300-350  K  and  thus  within  experi¬ 
mental  error  an  approximate  constant  value  can  be  assumed  as 
reasonable. 

6.1.  Application  of  the  model 

Recently  we  have  reported  an  improved  performance  DMFC 
operating  under  ambient  air  conditions  [19].  The  improved  per¬ 
formance  was  based  on  an  improved  MEA  fabrication  procedure 
using  a  carbon  paper  backing  layer.  We  have  applied  the  model 
discussed  in  this  paper  to  the  data  as  shown  in  Figs.  9  and  10. 
As  can  be  seen  there  is  good  agreement  between  the  model  and 
experiment. 

The  Tafel  slope  and  the  cell  resistance  (parameters  b  and  R 
in  the  proposed  empirical  equation)  were  kept  constant  and  the 
same  for  the  cell  tested  in  the  present  paper.  Although  the  perfor¬ 
mance  of  the  cell  with  the  modified  MEA  is  significantly  better, 
the  cathode  pressure  is  ambient  and  methanol  concentration  is 
higher  (1.0  M)  the  model  still  describes  the  cell  behaviour  quite 


298  308  318  328  338  348  358  368 


ke ff  =  k\0Cy[Ea  exp 


Cell  Temperaure/  K 


where  k\0 ,  a  and  A  are  constants  with  the  following  values: 


(10) 


Fig.  8.  Empirical  equation-based  predicted  values  for  the  kinetic  related  param¬ 
eter  a  as  a  function  of  cell  temperature  for  different  methanol  solution  concen¬ 
trations  (0:  0.125  M;  x:  0.25  M;  A:  0.5  M;  Q:  0.75  M). 
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Current  Density/  A  cm 2 

Fig.  9.  Demonstration  of  the  empirical  equation  to  an  oxygen  fed  DMFC  with 
different  ME  A  fabrication  technique.  The  cell  is  operated  at  ambient  pressure 
with  1.0 M  aqueous  methanol  solution  (cell  temperature:  ♦:  60  °C;  A:  90  °C). 


methanol  concentration  and  effective  mass  transport  coefficient, 
the  values  are  higher  for  the  new  cell  design.  For  example  at  60 
and  90  °C  the  values  are  0.9  and  0.45,  respectively  (for  0.75  M 
methanol),  with  the  original  cell  design  compared  to  2.5  and 
1.05  for  the  new  cell  design.  That  is,  effective  mass  transport 
coefficients  for  the  new  design  are  approximately  twice  those  of 
the  original  cell  design. 

The  values  of  C\  at  90  °C  are  similar  for  both  sets  of  data  at  the 
two  temperatures.  The  variations  in  the  values  of  C\  at  60  °C  are 
more  difficult  to  explain,  i.e.  approximately  0.1  for  the  original 
cell  design  with  air  and  now  0.085  for  air  and  0.058  for  air 
and  oxygen,  respectively.  However,  the  term  C\  is  essentially 
a  multiplier  on  the  mass  transport  limiting  component  of  the 
empirical  equation  and  thus  variations  in  the  values  are  indicative 
of  differences  in  cell  design. 

7.  Open  cell  voltage  modelling 

In  theory  the  open  circuit  potential  (or  rather  the  equilibrium 
potential)  for  the  DMFC  is  given  by  the  Nernst  equation: 


0.9 
0.8 
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>  0.6 
Q> 

|  0.5 
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>  0.4 
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U  0.3 
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Fig.  10.  Demonstration  of  the  proposed  empirical  equation  to  an  air  fed  DMFC 
with  different  MEA  fabrication  technique.  The  cell  is  operated  at  ambient  pres¬ 
sure  with  1.0 M  aqueous  methanol  solution  (cell  temperature:  ♦:  60  °C;  A: 
90  °C). 

accurately.  The  values  of  the  coefficients  used  to  fit  the  data  are 
given  in  Table  1 . 

Comparing  the  parameters  C\  and  C2  for  both  sets  of  experi¬ 
mental  data  we  see  that  there  are  some  significant  differences.  In 
the  case  of  C2  which  is  inversely  proportional  to  the  product  of 


E  _  r0  ,  (  <2ch3qh(};02  ,airPcath)  ^ 

°C  “  6F  n  1  ^ 


(11) 


where  Eoc  is  the  cell  voltage  (V),  the  standard  cell  voltage 
(V),  <2ch3oh  the  activity  of  methanol  in  solution,  yo2,air  the  mole 
fraction  of  oxygen  in  air  (— ),  pcath  the  cathodic  gas  pressure  (Pa), 
and  pco2  is  the  partial  pressure  carbon  dioxide  (Pa). 

However,  it  is  known  that  as  well  as  interfacial  phenom¬ 
ena  affecting  the  vale  of  Eq  a  major  factor  is  that  of  methanol 
crossover  which  depolarises  the  cathode.  This  causes  a  signif¬ 
icant  reduction  in  potential  from  that  predicted  by  the  Nernst 
equation.  As  a  model  of  the  cathode  of  the  DMFC  we  base  this 
on  a  mixed  potential  for  oxygen  reduction  and  methanol  oxida¬ 
tion. 

As  a  basis  we  adopt  Butler  Volmer  kinetics  for  methanol 
oxidation  and  oxygen  reduction  in  the  following  form,  which 
were  used  as  a  basis  for  the  empirical  non-equilibrium  model: 


.  (Cme)^ 

Jme  =  Jo  f  exp 
CME 


"aa^a" 

Pc  02 

PYD 

"-aaF?7a" 

RT 

CA 

PcOi  ref 

RT 

(12) 


J02  =  70c 


(po)N° 


(13) 


Table  1 

Model  parameters  for  the  improved  fuel  cell  performance 


02  (60  °C) 

02  (90  °C) 

Air  (60  °C) 

Air  (90  °C) 

£o(V) 

0.39 

0.395 

0.33 

0.34 

b  (V  dec-1) 

0.11 

0.124 

0.11 

0.124 

R  (£2  cm2) 

0.192 

0.0225 

0.192 

0.0225 

Ci(V) 

0.058 

0.12 

0.085 

0.126 

C2  (cm2  A-1) 

2.5 

1.05 

2.5 

1.05 

where  jja  =  E0c  -  £me  and  rk  =  E0c  -  £q2  and  pco2  refers 
to  pressure  of  carbon  dioxide  product  and  Cme  refers  to  the 
methanol  concentration  at  the  cathode  which  will  be  deter¬ 
mined  by  the  overall  rate  of  transport  of  methanol  through  the 
MEA. 

We  now  consider  that  the  total  current  at  the  cathode  is  the 
sum  of  the  methanol  oxidation  and  oxygen  reduction  currents: 

j  =  Jme  +  Jo2 


(14) 


K.  Scott  et  al.  /  Journal  of  Power  Sources  161  (2006)  885-892 


891 


Then  apply  equilibrium  such  that 

(Po)N°  ( ___  (  0CcFr]c 


7  =  0  =  70c 


ref 
P  o 


exp 


RT 


exp 


/  acFrja 
V  RT 


+70 


.  (CmeU 


z^-ref 

CME 


-exp 


otaFr]a 

Pc  o2 
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'-aaFrj  a' 

RT 

CA  JJ 

^C02ref 

RT 

(15) 


With  the  assumption  that  the  partial  pressure  of  carbon  dioxide 
will  be  small  we  simplify  this  to 


70c 


( poff 

ref 

P  o 


exp 
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.  (Cme)j 

-^o^^exP 

CME 


WqFQc 

RT 

RT 


exp 


OcPjh 
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(16) 


Re-arranging  and  taking  logs  gives 
ln(exp(-^ME^a)(exp(-^c^c)  -  exp(/3c^c))) 

_  !  7o/of  I  AT  1  f  6^Me\ 

where  generally,  ft  =  anF/RT. 

Hence 


(17) 


(-0ME??a)  +  ln(exp(-^c/7c)  -  exp(^c/7c)) 

,  jopro  |  tlj  i  f  Cme\ 

-  STpo?1" +Wm  Ic^J 


(18) 


which  gives 


(-^ME?7a)  -  PcVc  +  ln(l  -  exp(2/lc?7c)) 
,  jopro  |  tlj  i  f Cme\ 


(19) 


Under  conditions  when  the  cathode  is  significantly  polarised 
by  methanol  (which  occurs  in  practice)  then  from  the  definition 
of  overpotential,  noting  that  values  of  r]c  are  negative,  we  thus 
obtain: 

:  „  ref 

(Pme  +  Pc)E0c  =  PmeE^e  +  pcEf  -  In-  70/  0 


j0c(Po)N° 


VlME/ 


(20) 


In  fact,  several  researchers  have  observed  that  the  open  circuit 
voltage  decreases  with  increasing  methanol  concentration  and 
that  the  cathode  electrode  performance  was  significantly  lower 
at  higher  methanol  concentration  [4,18]. 

The  model  Eq.  (5)  states  that  the  open  circuit  potential  varies 
linearly  with  the  log  of  methanol  solution  concentration  accord¬ 
ing  to 


(21) 
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Fig.  11.  Experimentally  measured  and  model  predicted  (dotted  lines)  open  cell 
voltage  as  a  function  of  cell  temperature  (shown  in  figure)  for  a  cell  operated 
with  various  methanol  concentrations  supplied  at  a  rate  of  1.12  cm3  min-1  with 
air  fed  cathodes  pressurised  at  2  bar. 


where 


co  = 


PmE  0  ,  Pc  pO 

Pc  +  PmE  ME  Pc  +  PmE  °2 


_ 1 _ in  ( _ Ml _ ) 

Pc  +  Pme  y  70co (po)N°CmE  ) 


and 

l  =  -+L 

PmE  +  Pc 

which  vary  with  temperature. 

Fig.  11  shows  a  typical  fit  of  Eq.  (21)  to  experimental  data 
at  a  fixed  temperature.  Agreement  between  model  and  data  is 
quite  reasonable  with  the  coefficients  shown  on  Fig.  12.  The 
coefficient  X  varies  linearly  with  temperature  and  at  343  K  the 
value  is  —0.02755.  Hence  from  this  the  combined  value  of  Tafel 
slopes  as  defined  by  Z?ceii  is  0.63  which  is  comparable  to,  but 
smaller  than,  the  value  shown  in  Fig.  3  (assuming  Ame  =  1.0). 
There  are  a  number  of  factors,  which  would  contribute  to  this 
disagreement  in  derived  parameters  but  major  factors  will  be  the 
amount  of  experimental  data  and  assumptions  in  the  model. 
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Fig.  12.  Variation  of  the  fitting  coefficients  for  Eq.  (6)  as  a  function  of  cell 
temperature  for  the  data  used  in  Fig.  11. 
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Fig.  13.  Comparison  of  the  model  based  predicted  open  cell  voltage  with  actual  experimental  values  obtained  as  a  function  of  cell  temperature  and  for  various 
methanol  solution  concentrations  supplied  at  a  rate  of  1.12  cm3  min-1  with  air  fed  cathodes  pressurised  at  2  bar. 


The  coefficient  co  also  varies  linearly  with  temperature  (over 
the  range  considered)  as  predicted  by  the  model.  Overall  the 
variation  in  open  circuit  potential  as  predicted  by  the  model  is 
shown  in  Fig.  13. 

8.  Conclusions 

Experimental  data  from  a  small-scale  DMFC,  operated  over 
a  range  of  temperatures  and  moderately  low  methanol  solution 
concentrations,  has  been  used  to  validate  a  semi-empirical  model 
of  the  cell.  The  model  parameters  or  coefficients  have  been  deter¬ 
mined  by  the  best  statistical  fit  to  experimental  data  and  have 
been  used  to  determine  approximate  kinetic  and  mass  transport 
related  constants.  The  open  circuit  potential  of  the  cell  has  been 
interpreted  by  a  model  based  on  a  mixed  potential  associated 
with  oxygen  reduction  and  methanol  oxidation  at  the  cathode. 
The  semi-empirical  model  of  the  DMFC  has  been  used  to  model 
data  from  a  new  cell,  which  gives  significantly  higher  power 
densities  than  the  test  cell  used  to  validate  the  model. 

Acknowledgements 

The  authors  would  like  to  acknowledge  the  following  for 
support  of  this  research:  the  European  Commission  for  a  TMR 
Marie  Curie  B20  and  an  IHP  Marie  Curie  B30  research  training 
grant  to  Dr.  P.  Argyropoulos;  EPSRC  supported  Mr.  C.  Jackson 
through  a  PhD  studentship. 

References 

[1]  K.  Scott,  W.M.  Taama,  J.  Cruickshank,  Performance  of  a  direct  methanol 
fuel  cell,  J.  Appl.  Electrochem.  28  (3)  (1998)  289-297. 

[2]  K.  Scott,  W.M.  Taama,  J.  Cruickshank,  Performance  and  modelling  of  a 
direct  methanol  solid  polymer  electrolyte  fuel  cell,  J.  Power  Sources  65 
(1-2)  (1997)  159-171. 

[3]  S.F.  Baxter,  V.S.  Battaglia,  R.E.  White,  Methanol  fuel  cell:  anode,  J.  Elec¬ 
trochem.  Soc.  146  (2) (1999)  437-447. 

[4]  K.  Sundmacher,  K.  Scott,  Direct  methanol  polymer  elec¬ 
trolyte  fuel  cell:  analysis  of  charge  and  mass  transfer  in  the 


vapour-liquid-solid  system,  Chem.  Eng.  Sci.  54  (13-14)  (1999)  2927- 
2936. 

[5]  P.  Argyropoulos,  K.  Scott,  A.K.  Shukla,  C.  Jackson,  Investigation  of  the 
applicability  of  PEMFC  derived  empirical  equations  to  DMFCs.  A  new 
empirical  equation  approach  for  DMFCs,  Fuel  Cells-Fund.  Syst.  2  (2002) 
78. 

[6]  J.H.  Lee,  T.R.  Lalk,  A.J.  Appleby,  Modelling  electrochemical  performance 
in  large  scale  proton  exchange  membrane  fuel  cell  stacks,  J.  Power  Sources 
70 (1998)  258-268. 

[7]  J.H.  Lee,  T.R.  Lalk,  A  modelling  technique  for  fuel  cell  stack  systems,  in: 
Proceedings  of  the  ASME  Dynamics  Systems  and  Control  Division,  1996. 

[8]  G.  Squadrito,  G.  Maggio,  E.  Passalacqua,  F.  Lufrano,  A.  Patti,  An  empir¬ 
ical  equation  for  polymer  electrolyte  fuel  cell  (PEFC)  behaviour,  J.  Appl. 
Electrochem.  29  (1999)  1449-1455. 

[9]  J.  Kim,  S.  Lee,  S.  Srinivansan,  C.E.  Chamberlin,  Modeling  of  proton 
exchange  membrane  fuel  cell  performance  with  an  empirical  equation,  J. 
Electrochem.  Soc.  142  (8)  (1995)  2670-2674. 

[10]  P.  Argyropoulos,  K.  Scott,  A.K.  Shukla,  C.  Jackson,  A  semi  empirical 
model  for  the  direct  methanol  fuel  cell.  Part  1.  Model  development  and 
verification,  J.  Power  Sources  123  (2003)  190. 

[11]  P.  Argyropoulos,  K.  Scott,  W.M.  Taama,  Carbon  dioxide  evolution  pat¬ 
terns  in  operating  DMFC  cells,  Electrochim.  Acta  44  (1999)  3575- 
3584. 

[12]  P.  Argyropoulos,  K.  Scott,  W.M.  Taama,  Gas  evolution  and  power  per¬ 
formance  in  direct  methanol  fuel  cells,  J.  Appl.  Electrochem.  29  (1999) 
661-669. 

[13]  W.M.  Taama,  P.  Argyropoulos,  K.  Scott,  Parametric  studies  on  the  liquid 
feed  DMFC,  in:  Proceedings  of  the  1998  IChemE  Research  Event,  IChemE, 
Newcastle  upon  Tyne,  1998. 

[14]  K.  Scott,  W.M.  Taama,  P.  Argyropoulos,  Mass  transport  and  flow  visuali¬ 
sation  in  the  direct  methanol  fuel  cell,  in:  Proceedings  of  the  Electrochem 
99,  Portsmouth,  UK,  1999. 

[15]  K.  Scott,  W.M.  Taama,  P.  Argyropoulos,  Engineering  aspects  of  the  direct 
methanol  fuel  cell  system,  J.  Power  Sources  79  (1)  (1999)  43-59. 

[16]  K.  Scott,  W.M.  Taama,  P.  Argyropoulos,  Material  aspects  of  the  liquid 
feed  direct  methanol  fuel  cell,  J.  Appl.  Electrochem.  28  (12)  (1999)  1389— 
1397. 

[17]  K.  Scott,  P.  Argyropoulos,  Flow  and  mass  transport  in  direct  methanol  fuel 
cells,  Int.  J.  Transport  Phenom.  3  (2001)  257. 

[18]  M.K.  Ravikumar,  A.K.  Shukla,  Effect  of  methanol  crossover  in  a  liquid 
feed  polymer  electrolyte  direct  methanol  fuel  cell,  J.  Electrochem.  Soc. 
143  (8) (1996)  2601-2606. 

[19]  A.  Shukla,  K.  Scott,  C.  Jackson,  An  improved  performance  liquid  feed 
solid  polymer  electrolyte  direct  methanol  fuel  cell,  Electrochim.  Acta  111 
(2002)  43. 


